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Abstract 

The ATLAS collaboration published supersymmetry limits based on up to 
about 4.7 fb _1 data collected over the year 2011 from LHC runs at 7 TeV. 
These were mainly interpreted within restricted, particular or simplified mod- 
els for supersymmetry breaking schemes or scenarios. The pMSSM is an al- 
ternative and more generic supersymmetry framework which captures broader 
phenomenological features. Searching for more generic conclusions from the 
supersymmetry limits interpretation, we update a Bayesian global fit of the 
pMSSM to pre-LHC data using the LHC-7 limits. The posterior distributions 
show the most up to date features, revealing allowed versus excluded regions 
in sparticle mass planes within the MSSM. 



1. Introduction 

The discovery of Higgs boson-like state around 126.5 GeV 1 1 ] or 125 GeV (2J at the large hadron collider 
(LHC) is an excellent accomplishment. However this only marks the beginning of exciting moments for 
particle physics endeavour in establishing the mechanism of electroweak symmetry breaking and for 
shedding light on new physics beyond the Standard Model (BSM.) During the year 2011 run of the 
LHC machine both ATLAS and CMS detectors have recorded up to about 5 fb _1 of data. The collab- 
orations have conducted many analyses on the data, searching for, amongst other new physics models, 
supersymmetry (SUSY) by looking for final states containing jets and large missing transverse energy 
(MET) that could indicate the production of squarks and gluinos in the collider |3||4j. There is no sign 
for SUSY observed as of the time of writing this article according to public results. As such findings of 
the experiments were presented in the form of model-independent non-SM cross section limits and inter- 
pretations showing exclusion regions within specific models of SUSY such as the constrained version of 
the R-parity conserving minimal supersymmetric standard model (MSSM) called CMSSM/mSUGRA, 
particular SUSY breaking schemes (for a review see e.g. [5]) or simplified models [6] of SUSY scenarios. 



1.1 Aim of the article 

In this article we will assess the impact of the LHC-7 SUSY results on the R-parity conserving phe- 
nomenological MSSM (pMSSM). In particular we are going to use the ATLAS SUSY limits reported in 
Refs. (7HT7J to discriminate between allowed and excluded regions in the pMSSM sparticle mass planes. 
To date several research groups have looked into analysing the impact of LHC-7 data on various SUSY 
models. For a non exhaustive instances see Refs. 1 18 -34). Out of these, the interpretation, within the 



pMSSM, for the CMS collaboration's SUSY limits based on 1 fb x data presented in ||22p has a particular 
relevance. Our analyses updates the impact of LHC results on the pMSSM, here by using SUSY limits 



from the ATLAS collaboration with up to about 4.7 fb x data. Following the approach in Ref. [22], the 



data d from the various experiments used for our analysis is decomposed into two independent parts, 

d = d pre _LHC © d-LHC (1) 

where d pre -LHC represents the indirect pre-LHC collider and cold dark matter relic density constraints 
summarised in Table [T] and d^uc, the LHC-7 SUSY limits shown in Table [2] The pMSSM posterior 
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distributions from the Bayesian global fit [36] to pre-LHC data is now used as the prior, n(6), for updating 
the pMSSM posterior sample with the LHC SUSY limits. Using Bayes theorem, weighing the prior with 
the likelihood over the LHC data L{diHc\G) gives an updated, post-LHC, posterior distribution 

p{0\dLHc)~L{d LH c\0)*{p) (2) 

valid up to a normalisation factor. Analysing this will reveal information about the impact of the LHC-7 
data on the pMSSM parameter space. 

This article is structured as follows. In remaining part of this Section we give a brief recapitulation 
of the Bayesian global fit of the pMSSM to the pre-LHC data allowing us to set the context for explaining 
our analysis. In Section [2] we present the methodology employed in simulating SUSY event at the 
LHC and the computation of the pMSSM predictions for the cross section within acceptance, the main 
observable to compare with the upper limits from ATLAS. The impact of the experimental limits on the 
pMSSM is presented in Section [3]J The last Section is preserved for discussing our conclusions and 
outlook. 



1.2 The pre-LHC pMSSM fit review 

The Bayesian global fit of the pMSSM to pre-LHC data were performed in [35][36| . The posterior 
samples reveal SUSY spectra with various characteristics satisfying different phenomenological scenar- 
ios |37] - [39| which are mainly difficult to capture within the classic constrained benchmark models. For 
the pMSSM fit, the parametrisation is completely decoupled from the details of the physics responsible 
for SUSY breaking. Requiring compatibility with observations about CP violation and flavour changing 
neutral current processes, only real soft SUSY breaking terms are considered, with all off-diagonal ele- 
ments in the sfermion mass terms and trilinear couplings set to zero, and the first-and second-generation 
soft terms equalised; leading to a set of 20 parameters: 

= {M 1)2 , 3 ; rn^ DLy mfj^- A tjb , T ,, =e , tan/3}, (3) 

where M\, M2 and M3 are the gaugino mass parameters; and mj are the sfermion mass parameters. 
A tj b ,r,[i=e represent the trilinear scalar couplings while the Higgs-sector parameters are specified by the 
two Higgs doublet masses mjy , Tn 2 H , the ratio of the vacuum expectation values tan fj = (H2) / (Hi) 
and the sign of the Higgs doublets mixing parameter, sign(fj,). The pre-LHC data, d pre -LHC = §_ = 
{fj,i, (Xi} where % = 1,2, . . ., represents the experimental central values and errors for the electroweak 
physics observables, B-physics observables and the cold dark matter relic density, summarised in Table[T] 

O = {m w , sin 2 9 l e J f , T z , <5a M , 2$, A%, A 1 = A e , H° c , A b f b , A b > c , (4) 
BR{B^X sl ), BR(B S — > fi + fjT), A _, R B r(b u ^tv), «aa/ Ss , 

DM h }. 

The pre-LHC posterior distribution from the Bayesian global fit which is now considered as a prior 
distribution for the update with the LHC data is given by 

tt(0) = tt'(0) J] (2^ 2 )~ 1/2 exp [-(0« - inf/2al] (5) 
i=l 

where ir'(Q) is the prior probability density for the Bayesian global fits to d pre _LHc which can be flat 
over the individual parameters Oi (for flat prior fits) or flat over the logarithm of the parameters, log 9 
(for log prior fits). Here we do not aim at checking the constraining strength of the LHC data over the 
pMSSM parameters. Thus no prior dependence analysis is discussed and we work with only the log 
prior posterior-samples of the pre-LHC global fits. In the next Section we describe the data, i.e. the 
extra-SM cross sections within acceptance, the simulation of SUSY events at the LHC, and the ATLAS- 
like analyses of the events. These are used in constructing the likelihood, L{dLHc\&) required in addition 
to the prior, tt(6), for completing the required variables in the target eq.Q. 
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Observable Constraint 


Observable Constraint 


m w [GeV] 80.399 ± 0.027 
T z [GeV] 2.4952 ± 0.0025 
sin 2 6 l e J f 0.2324 ± 0.0012 
Sap (30.2 ± 9.0) x 10 10 
Itf 20.767 ± 0.025 
E? b 0.21629 ±0.00066 
R° c 0.1721 ±0.0030 
Apg 0.0992 ± 0.0016 
A^ B 0.0707 ± 0.035 


A 1 = A e 0.1513 ±0.0021 
A b 0.923 ±0.020 
A c 0.670 ±0.027 
Br(B -»• X s j) (3.55 ± 0.42) x 10 4 
Br(B s ->■ < 5.8 x 10~ 8 (see caption.) 
i?AM Bs 0.85 ±0.11 
RbhbJtu) 1.26 ±0.41 
A - 0.0375 ± 0.0289 
ttcDhih 2 0.11 ±0.02 



Table 1: Summary for the central values and errors for the electroweak physics observables, B-physics observables and cold 
dark matter relic density constraints. The posterior distribution from the pre-LHC fit for Br(B s — > fi + fi~) is centred around 
2.8 x 10~ 9 so most of the points are in agreement with the more recent bound 4.5 x 10 -9 |40|. 



2. Analysis 

Our analyses is centred around computing the likelihood for dmc given the pMSSM parameters, f?j. 
In this Section we describe the ATLAS data and the simulation of SUSY events. The latter allows for 
computing the extra-SM, here SUSY, cross sections within acceptances, a acc , over various cuts on the 
collider final state characteristics. The degree of agreement or deviation between the pMSSM predictions 
for a acc and the experimental values is used to quantify the plausibility of obtaining the data from the 
model parameters, L(d,LHc\Q)- 



2.1 The data, d LHC 

The LHC data we use are the 95% C.L. upper limits on the non-SM cross section within acceptance. 
For each of the ATLAS analyses (7- 17 1, there are various signal regions defined by specific set of cuts 
and events selection criteria. The results are based on various, namely 35 pb _1 up to 4.7 fb _1 , data set 
recorded by the ATLAS detector at 7 TeV centre of mass energy run of the LHC in the year 2011. The 
analyses were designed to capture different possible manifestations of SUSY after the proton-proton 
collisions at the LHC. 

SUSY production at the collider would be dominated by large direct production of squark and 
gluino pairs (g g, g q, or q q) that would decay (q — > qxi and g — >■ qq Xi) to the weakly interacting 
neutralino lightest supersymmetric particle (LSP), x?> which escapes the detector unseen in the form 
of missing transverse energy, MET. The different groups of search channels, which we briefly describe 
here, are all MET-based. The first is the search for squarks and gluinos that lead to final states containing 
high-pT jets, MET and no leptons (electrons or muons) as in Refs. (7] [TT]jT4| . The strategy for this group 
of searches is optimised for maximal discovery reach in the rrig-mq plane. This group of search channels 
could be specialised to the case of having heavy flavour jets. Doing this would capture the scenario 
where the sbottoms (b\) and stops {t\) are lighter than other squarks such that direct or gluino-mediated 
production (g — > bb or g — > ti) is the dominant SUSY production mode in the collider as considered in 
Refs. QDJ[TS). 

Requiring final states containing one (or more) electron or muon in addition to jets and MET would 
capture scenarios where gluinos cascade decay products include a charginos, which subsequently 
decays into final states containing high-py leptons as considered in Refs. ||9jfT5l. Further, in scenarios, 
such as Natural SUSY |4T[ , where first and second generation sfermion masses are larger than few TeVs, 
the direct production of weak gauginos may be the dominant SUSY processes. When both gauginos 
decay leptonically, a distinctive signature with no jets, three leptons and significant MET, as considered 
in Ref. (TTJ, can be observed. 
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Channels + MET 


Signal rej 


jions, &bsm upper limits 


Luminosity 


jets + O-lep 


1.3 


0.35 


1.1 


0.11 


- 


- 


- 


- 


- 


- 


- 


35 pb- 1 [ 


71 

[11] 


jets + O-lep 


22 


25 


429 


27 


17 


- 


- 


- 


- 


- 


- 


1.04 ft)- 1 


>6 jets 


194 


8.4 


12.2 


4.5 


- 


- 


- 


- 


- 


- 


- 


1.34 ft- 1 


[12] 

*P 


Ojet + 2-lep 


0.22 


0.09 


0.21 


0.07 


0.07 


0.07 


- 


- 


- 


- 


- 


35 pb- 1 [ 


jets + 1-lep 


41 


53 




















165 pb" 1 !^ 


T 


b-jet + O-lep 


288 


61 


78 


17 
















830 pb- 1 [1( 


>(2-6)jets + O-lep 


0.62 


5.3 


6.2 


0.65 


3.5 


3.7 


12 


2.2 


2.6 


2.5 


18 


4.7 ft- 1 [ 


13 


i 


>(6-9)jets + O-lep 


14 


4.2 


1.2 


9.8 


3.2 


0.81 












4.7 ft- 1 [ 


14 


i 


>(2-4)jets + 1-lep 


1.3 


1.5 


3.7 


















4.7 ft- 1 [ 


15 


i 


>(l-2)b-jets + O-lep 


283 


65 


15.4 


61 


14.4 


4.3 


22.2 


8.5 








2.05 ft- 1 [16] 


3-lep 


3.5 


1.5 




















2.06 ft- 1 [17] 



Table 2: The ATLAS 95% C.L. upper limits on the extra-SM cross sections within acceptance for the various signal regions 
described in the text. The limits on each search channel row are ordered; with the first representing the first corresponding 
name of the signal region described in the corresponding experimental paper. The unit for each cross section is the inverse of 
the corresponding luminosity. 



In all, we considered 55 SUSY signal regions. For each region the number of events that pass 
the selected criteria and also the expected Standard Model (background) events were reported Q In 
addition the upper bound on the cross sections for non-SM interactions were also given. These allow for 
comparisons with any BSM predictions to determine whether the new physics model is allowed or ruled 
out at the 95% confidence level. This is the approach we have chosen. The LHC data, (Ilhc, for our 
analysis is represented by the set summarised in Table [2] The limits are used to constrain the pMSSM 
o" acc predictions from the simulation of SUSY production at the LHC. 



2.2 The pMSSM predictions for the cross section within acceptance, of cc 

In order to compute the predictions for af cc within acceptances over the cuts and selections criteria that 
defines the various ATLAS SUSY signal regions, we simulate the generation of SUSY events at 7 TeV 
LHC using Monte Carlo events generator and then analyse the collider final states in a similar fashion to 
the ATLAS procedures. 



2.2.1 SUSY events simulation 

We use Herwig++ [42] 43 1 to simulate four sets of sparticle production processes for each point in the 
pMSSM sample, namely: a) squark-squark and squark-gluino production, b) the production of an elec- 
troweak gaugino in association with a squark or gluino and c) the production of slepton and electroweak 
gaugino pairs. Each of the pMSSM posterior sample point in the SUSY Les Houches Accord (SLHA) 



file format produced by Sof tSUSY (44) is passed to Herwig++ [|42j|43[ for generating 1000 SUSY 
events. Through out the analysis we use the SUSY production cross section from the event generator 
calculated at leading order in perturbative QCD0 



2.2.2 Analysing the simulated SUSY events using Rivet 

The analysis of the generated SUSY events are done at the particle level using the "Robust Independent 
Validation of Experiment and Theory", Rivet, Monte Carlo validation framework (45j. We use this 

'No detector simulation is attempted as that would make no significant effect on our findings or conclusions. 
2 About half of the posterior samples we consider have a negative gluino mass, allowed by the SLHA accord, that crashes 
the NLO calculator we have access to. As such the NLO correction is dropped out. 
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No. 


Constraint 


Log prior survive 


1. 


m h = 122 - 128 GeV 


36.08% 


2. 


m h = 125 - 126.5 GeV 


9.17% 


3. 


SUSY a acc limits 


57.47% 


4. 


m h (l) and SUSY a acc limits 


14.53% 


5. 


1.67-CT i?ZZ;77 


13.30% 


6. 


mh(l) & 1.67-cr Rzz m 


4.67% 


7. 


1.61 -a Rzz-rn & SUSY <j acc limits 


6.11% 


8. 


m h , 1.67-cr i? z z; 77 & SUSY a acc limits 


2.19% 



Table 3: Summary of the relative number of surviving posterior points, from the pre-LHC Bayesian global fits of the pMSSM, 
and after imposing the Higgs discovery data and SUSY limits. 



to analyse each and every SUSY event generated by the Monte Carlo collider simulator, without the 
need for detector simulations, and the publicly available Rivet analyses for the ATLAS SUSY searches 
in Refs. [|7j-[T7]l. The Rivet analyses are plugged-in to Herwig++ for computing the acceptances, Ai = 
N cu ts/NtotaU after applying the various cuts on the kinematic variables of the collider final states. Here 
N cuts is the number of events that pass the experimental cuts and N to tai = 1000 is the total number 
of generated SUSY events. Rivet acts per-event wise on the events produced by Herwig++. The jet 
identification is done using fast jet p6j . The cross section within acceptance is computed as 

acc c * SUSY,LO ( ~ 

a i = eAi a i (6) 

where we consider the efficiency, e = 1 (since no detector simulatiorj^]); i = 1, 2, . . . , 55 runs over the 
55 different signal regions, and a SUSY > LO j s the total LO SUSY production cross section. 



2.3 The likelihood, L(d LHC \9) 

The likelihood is a simple step function that equals to unity if the ATLAS limits are satisfied else zero if 
excluded. SUSY points with predicted cross sections smaller (greater) than the ATLAS non-SM limits 
are then allowed (excluded) at 95% confidence level according to 

55 

L(d LHC \0) = Yl£i; £i = 

i=l 



arc ^ acc,max *\ 

T >°i \ 

ncr ^ acc,max ( • 

•acc< f 



3. The SUSY limits' impact on the pMSSM 

The relative number of the pMSSM points that survive the SUSY limits and in various combinations 
with the Higgs-sector data (in the di -photon channels as applied in Ref. HTJ and other ATLAS, CMS, 
LEP and Tevatron Higgs-sector constraints implemented in the HEP packages HiggsBounds (48) and 
FeynHiggs |49j) are summarised in Table [3] The posterior probability distribution for the sparticle 
masses derived from the post-LHC distribution eq.Q are shown in FigjT] As can be seen from the 
mentioned plots, the ATLAS SUSY limits are most sensitive in constraining the gluino mass whose 
central value is now shifted^ from around 2 TeV to 3 TeV followed by the squark masses. The effect 
of the SUSY limits on the gluino mass seems to be very different from what happen for the case of 
constrained models such as CMSSM/mSUGRA (addressing this is outside the purpose of this article). 
The limits are relatively less sensitive in constraining sleptons, electroweak gauginos and the lighter 

Performing the analysis with fast detector simulation will not have a significant effect on our results or conclusions. For 
instance, a ±2 GeV accuracy gain in the MET or jet pr will have little or no affect on our results. 

4 The SUSY search channels that involves final states with high number of jet multiplicities severely constrains the gluino 
mass to higher values relative to the pre-LHC mass distribution. 
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sbottom or stop quarks. The resultant effect of applying the limits coming from the 55 different ATLAS' 
SUSY signal regions shown in Table [2] on the pMSSM can be summarised in plots showing allowed 
versus ruled out regions in sparticle mass planes. For gluino-stop and gluino-sbottom mass planes the 



plots are shown in Fig. 



Region where max{- 



acc.raax 



> 1.0 are excluded at 95% confidence level by 



the combined ATLAS fimits. The plots shows that unlike the gluino mass which is now constrained to 
be heavy ~ 3 TeV, the third generation squarks can be much lighter. 



4. Conclusions and outlook 

We have computed the effect of SUSY searches results from ATLAS collaboration with up to 4.7 fb _1 of 
7 TeV LHC data sets taken during the year 2011. This is done by simulating SUSY events at the LHC 
with Herwig++ (42j and the experimental analyses of the data with Rivet [45]. The particle-level 
per-event Rivet analysis is used for calculating the pMSSM predictions for the cross section within 
acceptances of the analyses cuts and events selection criteria. Comparing with the experimental 95% 
C.L. upper limits (shown in Table [2]) rule out about 40% of the initial pMSSM sample from a pre-LHC 
global fit to data. The SUSY limits are most sensitive in constraining the gluino mass whose distribution 
is now centred around 3 TeV (shifted from 2 TeV.) There is also a shift, but less significant compared 
to gluino mass case, in preference for heavier lst/2nd generation and stop masses as can be seen in the 
posterior distributions in Fig. [T] The excluded versus allowed regions by the ATLAS' extra-SM cross 
section upper limits on the gluino-stop and gluino-sbottom mass planes are shown in Fig. [2] 

Combining the SUSY and Higgs boson discovery data as described in Ref. [47] further constrain 
the pre-LHC posterior samples as summarised in Table [3] with only a single point (out of the initially 
about 40000 pre-LHC posterior sample points) surviving all of the following requirements: 



SUSY a acc upper limits, 

m h = 125 - 126.5 GeV, and 
Rzz- m = — = 0.56 ±0.25. 



Here /jlx = iFTj^ffr^ Br{h^>x] S M wnere = 77 or ZZ. The spectrumis characterised with a quasi- 
degenerate sbottom and LSP and heavy gluino and lst/2nd generation squarks. It is of the difncult-to-see 
at the LHC type discussed in Refs. [38,50]. It is worth mentioning that the data from the Higgs sector 



(4.67% model points survived) is far more constraining compared to the SUSY limits (57.47% model 
points survived). There is however some complementarity between the two set of constrained since 
applying the SUSY limits on the Higgs-data surviving models points brings down the surviving number 
to 2.19%. 

Our results go beyond the 1 fb _1 analysis done in Ref. fl22~| from various perspectives. First, our 
pre-LHC prior construction takes into account the constraint on the neutralino LSP relic density. We use 
the more stringent ATLAS SUSY limits from up to 4.7 fb _1 of data which include channels that better 
constrain gluino production with subsequent decay chains with several jets in the final states. 

The results and conclusion obtained from our analyses are conservative. This is because: the 
SUSY signal simulations were done at leading order (LO) in perturbative QCD as opposed to the next- 
to-leading order (NLO) results reported by ATLAS |^] Since the NLO cross sections are generally greater 
compared to the LO values, the exclusions here are more conservative, we will be allowing points that 
otherwise will be ruled out as is the case for the SPS4 bench mark point in Tab. [4] where the the LO 
and NLO cross sections within acceptance for snow mass points and slopes benchmark points can be 



5 About half of the volume of the pMSSM sample we employ have a negative gluino mass. This feature is the standard 
format allowed by the SLHA conventions. However such negative gluino SUSY points are not possible for the NLO cross 
section calculator that we have access to. 
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1 2 3 





1 2 3 






Fig. 1: The plots compare the log prior pMSSM sparticle masses' marginalised 1 -dimensional pre-LHC posterior distributions 
(dashed-red curves) and the surviving parameter regions after imposing only the SUSY limit (black curves) and both rrih = 
122.0 — 128.0 GeV and SUSY limits together (blue curves). All the masses are in TeV units. The vertical axes represent the 
relative probability weights of the model points. 
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M ffiuino [GeV] M gluiru > [GeV] 

Fig. 2: The plots show 95% confidence level exclusion contours for gluino-stop and gluino-sbottom mass planes derived from 
the combined set of ATLAS limits. The colour scales are proportional to the expected number of signal events normalised to the 
combined exclusion limit. The contour max{ —^rtojj ) = 1-0 determines the exclusion boundaries within the sparticle masses 
plane. Region with colour code greater than unity are excluded at 95% confidence level. 



compared. There is an approximate agreement, within the expected ~ 50% accuracies, with the corre- 
sponding values obtained by previous computations (with NLO corrections) (27J. 

For the analysis here only the posterior samples from a pre-LHC global fit to data with a loga- 
rithmic prior distribution over the 20 pMSSM parameters were considered. No analysis is done with 
the flat prior sample because it know that the pre-LHC fits were prior-dependent. However, it will be 
interesting [52| to estimate the strength of the LHC data by checking whether it allows for prior inde- 
pendent results which is necessarily needed for making conclusions regarding the predictive power of 
the pMSSM. This seems possible given the apparent interplay between Higgs boson decay rate in the 
di -photon decay channels which would require light sparticles^Jand the absence of SUSY signal to date 
at the LHC. 
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Benchmark point 


cr/pb 


status 




A 


B 


c 


D 


ATLAS 35pb _i 


ATLAS limits 
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0.11 
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1.347 
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1.172 


0.299 


A,B,C,D 
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0.933 
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0.418 


A,B,C,D 
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0.057 


0.062 


0.041 


allowed 
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0.089 


0.098 


0.067 


allowed 
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0.499 
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0.425 


0.169 


D 
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